Resistivity saturation revisited: results from a dynamical mean field theory 
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We use the dynamical mean field method to study the high-temperature resistivity of electrons 
strongly coupled to phonons. The results reproduce the qualitative behavior of the temperature 
and disorder dependence of the resistivity of 'A- 15' materials which is normally described in terms 
of saturation, but imply that the resistivity does not saturate. Rather, a change in temperature 
dependence occurs when the scattering becomes strong enough to cause a breakdown of the Migdal 
approximation. PACS: 72.10-d,72.10.Di,72.80.Ga 



The conventional |jj theory of the temperature depen- 
dent resistivity of metals predicts that at temperatures 
greater than the phonon Debye frequency ojd the resis- 
tivity, p, behaves as 



A p T - 



B„ 



(1) 



with the temperature coefficient A p proportional to the 
electron-phonon coupling strength and residual resistiv- 
ity Bp proportional to the impurity concentration. The 
theory accounts well for the properties of many metals 
but fails in a number of important cases. The first to at- 
tract widespread attention were the 'A-15' materials such 
as ViSi and Nb^Ge. In these compuinds the electron- 
phonon coupling is unusually strong and although the 
high temperature p(T) may be crudely fit to Eq. 1, A p 
is much smaller than expected and, crucially, B p is large 
even in pure samples Further, the effects of extra 
disorder (induced e.g. by radiation damage) are not ad- 
ditive: the extra disorder-induced resisitivity is temper- 
ature dependent, being larger at lower T and smaller at 
higher T; indeed for sufficiently high disorder levels p(T) 
becomes flat and then dp/dT becomes negative ||. This 
complex of behaviors was named 'resistivity saturation' 
by Fisk and Webb M . Interest in this long-standing prob- 
lem has been increased by observations that in a number 
of materials in which electron-electron interactions are 
believed to be important, including the high-T c cuprates 
|| , the A 3 C 60 materials @, and SVRu0 3 §, p(T) in- 
creases rapidly with increasing temperature and exhibits 
no signs of saturation even in the sense defined above, 
although the resistivity becomes much larger than that 
of A-15s. A recent paper argued that this absence of 
saturation was itself anomalous and interesting. 

Ref Q has led to a large but still inconclusive liter- 
ature on the cause of the phenomenon p|-pd|. Eq |l| 
is the mathematical statement that p is proportional 
to lattice disorder, i.e. that electrons are scattered by 
displacements, due either to thermal fluctuations or de- 
fects in the crystal, of ions (labelled by i) from their 
ideal crystallographic positions. Its derivation is based 



on three assumptions: (1) at temperatures greater than 
ujd, the thermal fluctuations of the ions may be modelled 
in terms of classical harmonic oscillators with spring con- 
stant k; equipartition then implies (xf) ~ T/k plus a T- 
independent term arising from defects, (2) the electron- 
ion interaction may be treated by second-order pertur- 
bation theory so that if the electron-ion coupling is g the 
electron scattering rate is g 2 (xfS and (3) usual Boltz- 
mann transport theory may be used to relate the resis- 
tivity to the scattering rate. Combining these consider- 
ations yields Eq |l| with A ~ g 2 /k and B proportional to 
the impurity concentration. 

Deviations from Eq [l] must come from the breakdown 
of one of the assumptions. Many authors have focussed 
on (3), however attempts to generate a systematic se- 
ries of 'multiple scattering' corrections to the Boltzmann 
equation have produced a multiplicity of terms whose in- 
terpretation and range of validity have not been clear 
(see, e.g. ||). However, examination of assumption (3) 
led Allen and Chakraborty to an important observation 
concerning the magnitude of the resistivity In the 
A-15 materials, the carriers are d-electrons and although 
the total carrier density is 4 d-electrons per metal ion, 
i.e. n = 12 per formula unit, the band structure is such 
that only one d-band crosses the fermi surface implying 
n = 1 per formula unit. The remaining d electrons are 
contained in filled bands which are very close to the fermi 
energy and to empty bands. A moderately large scatter- 
ing rate will mix the full and empty bands with the con- 
duction band, leading to a drastic increase in the effective 
number of carriers. Ref || concluded that while the effect 
does not fully explain the observed 'saturation' behavior, 
it does permit the magnitude of the high temperature re- 
sistivity to vary dramatically between different systems. 
We note here that the ruthenates and cuprates do not 
have many filled and empty d-bands near the fermi sur- 
face, and so may have a lower carrier density and thus 
higher resistivity than do the A-15s. 

Belitz and Schirmacher Q used a 'memory function' 
formalism to study a model of electrons interacting with 
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phonons via both density and stress couplings. In their 
model p increased indefinitely with T due to an interplay 
between scattering caused by the density coupling and 
hopping induced by the stress coupling. Our results sug- 
gest that the stress coupling and associated hopping are 
not necessary to reproduce the phenomenon. 



Yu and Anderson |11| focussed instead on assumptions 
(1) and (2), arguing that a strong electron-phonon cou- 
pling could lead to a spontaneously generated double- 
well potential for the phonon, rendering the second-order 
perturbation theory treatment invalid and requiring a so- 
phisticated quantum treatment. Our results support the 
idea that a strong electron-phonon coupling is the essence 
of the problem but imply that neither double-well forma- 
tion nor quantum effects are essential. 

In this Letter we calculate the resistivity of a model 
of electrons coupled with arbitrary strength to phonons. 
Because we seek to understand the fundamental mecha- 
nism underlying the saturation phenomenon we consider 
the simplest possible model: spinless electrons coupled to 
dispersionless phonons and to static disorder. Because 
saturation is a high-T phenomenon we further special- 
ize to classical phonons (cjn — > 0). The Hamiltonian 
H = H el + H ph + H e i- P h + H disorder with electronic part 
H e i — — Yip e pdpd p describing motion of electrons in a 
band with dispersion e p phonon part H p h = J2i \^ x \^ an 
electron-phonon coupling H e i- p h — gYi x i(d^ di — n), 
and an disorder part Hdisorder — Yi w i(df di — n). The 
form for Hdisorder corresponds to electrons scattered by 
random point defects which may be thought of as frozen- 
in lattice distortions of amplitude wi/g. The mean den- 
sity of electrons is n so x = is the equilibrium phonon 
state for a uniform distribution of electrons. We assume 
the random site energies Wi associated with the disorder 

are distributed according to Pdi S {w) = 
cxp(— w 2 / rj 2 ) / 2-Krj 2 ) . 

We define the parameter t to be one quarter of the full 
electron bandwidth; we measure all energies and temper- 
atures in units of t. Rescaling the phonon coordinate 
x —> x/Vk shows that the dimensionless parameter de- 
scribing the electron-phonon coupling is A = g 2 /(kt) (this 
A is a factor of it larger than the conventional MacMil- 
lan parameter). At low T the model may be solved by 
the usual Migdal pertubation theory, which is an expan- 
sion in Xy/ max(u>£>, T/t jL2j. To treat arbitrary coupling 
strengths we employ the dynamical mean field approxi- 
mation which becomes exact in the limit of spatial 
dimensionality d — ► oo and which comparisons to other 
techniques and to experiment have shown to be quan- 
titatively accurate in d — 3. The approximation may 
be derived from the assumption that the electron Green 
function G is 

G(c p , iuj n ) = [iuj n - E(iu n ) - e p + 1 . (2) 

where the self-energy £ is taken to be p- independent. 



This ansatz implies that all interaction effects are 
derivable from the local (p-integrated) Green func- 
tion, Gioc given by G ioc (iw„) = J G(e p ,iu) n ) = 

f -■ < ^p T> ( e p) _ q j s itself given in terms of a mean 

field function a(u>) determined by an equation which de- 
pends on the density of states T>\ we choose the semicir- 
cular form T> sernl 



a(uj) 



;(e p ) = ^At 2 - e 2 p /(2n t 2 ), so that 

Pph (x, w)P d is(w) 



dxdw- 



a(uj) + gx + w 
with the phonon probability distribution 
1 



(3) 



Pphonon — ^ €Xp 
^loc 



X 

2T 



+ / doj In [a(u>) + gx + w)] 



(4) 



with Z ioc = J dxdwP p honon(x,w)Pdis{w). G/ oc (w) = 

a a" ( fy c and £(w) = a(cu) - G^(u). This treatment 
of static disorder is equivalent to the familiar coherent 
potential 

We solve the equations numerically; computational de- 
tails are given in ref |Q. The conductivity may be cal- 
culated from |l3|] 

a ~ J da; J de p V(e p )[Im G(e p7 uj] 2 cosh~ 2 -^ (5) 

(Ref |D| takes the current operator to be independent of 
e p ; this does not affect results in any important way). 
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Fig 1 Resistivity p(T) vs temperature T (in units of 
t=bandwidth/4) at various couplings A for density n=l/4. Inset: 
Same couplings, density n=l/2. 

Fig. 1 shows the calculated resistivity p vs. temper- 
ature T at various couplings (assuming r) — i.e. no 
disorder) for density n=l/4; the inset shows the same 
calculation for n=l/2. Clearly the two fillings display 
the same qualitative behavior: at weak coupling, the re- 
sistivity has the usual linear T-dependence with a zero 



2 



offset and a slope proportional to A. (Quantum effects 
neglected here would cause the resistivity to drop dra- 
matically once T is reduced below u>d but do not affect 
fL2| the behavior at T > u>d; our results are only mean- 
ingful for T > lod)- For the narrow bands typical of 
A-15 materials, T = O.f corresponds to a temparture of 
order 300K, and only the results for T > 0.1 should be 
regarded as physically meaningful. At intermediate cou- 
pling, the high-temperature resistivity displays the es- 
sential features of the 'saturation' behavior found in the 
data, namely an apparently linear T dependence with 
non-zero offset and slope rather weakly dependent on A. 
At very strong coupling, the high-T behavior is not much 
changed (except that the offset becomes larger relative to 
the slope) but at low T the system becomes insulating 
(dp/dT < 0) below a A- dependent characteristic tem- 
perature Tg a p because a gap of size T gap opens at low T 
in the electron spectral function. Our interest is in the 
high-T behavior; the low-T insulating regime is exten- 
sively discussed in Ref JlJ]. 

We now determine the origin of the saturation behav- 
ior. The first issue is the relation between the resistivity 
and the electron self-energy. Fig 2 plots the calculated 
resistivity versus the imaginary part of the electron self- 
energy evaluated at the fermi surface, = 0)). 
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Fig 2 Resistivity p plotted against fermi surface scattering rate 
S (given in units of t=bandwidth/4) for various couplings, with 
temperature as an implicit parameter. The points which break 
away from the universal curve correspond to the low T insulating 
regime. 

One sees that except in the insulating low-T strong 
coupling gapped regime, p is a universal function of 
S"(0). At weak coupling, p ~ S"(0) as expected; at 
stronger couplings p increases faster than E"(0). Thus 
assumption (3) is not the issue. 

We next consider phonon anharmonicity. Fig 3 shows 
the mean square displacement of the oscillator coordinate 
as a function of temperature for the coupling strengths 



used in Fig. f .At weak coupling one sees the classical be- 
havior (a; 2 ) ~ T with zero intercept. As the coupling is 
increased the high-T behavior acquires the form (x 2 ) ~ 
AphononT + B phonon . However, a comparison of Figs 1 
and 3 shows the ratio B p ^ OTlorl / 1 A p h onon is much less than 
the ratio J5 p /A p ] in other words, the resistivity curves 
are much flatter than the lattice displacement curves. 
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Fig 3 Mean square lattice displacement < x 2 > vs. temperature 
T in units of t=bandwidth/4 at various couplings. 

Further, the calculated phonon probability distribu- 
tion P(x) (not shown) reveals that the double-well form 
discussed by Yu and Anderson occurs only in the low T, 
large g 'insulating' regime. Thus anharmonicity, while 
quantitatively important, is not the fundamental cause 
of the phenomenon. 

The remaining issue is assumption (2): the relation of 
the self-energy to the scattering mechanism. The weak- 
coupling (Migdal) result for our model is S (u> = 0) = 
g 2 < x 2 > (1 - (p/2t) 2 ) 1 / 2 . To see how this relation- 
ship evolves we plot in Fig 4 S (id 
function of temperature at different couplings 



0)/g 2 < x 2 > as a 
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Fig 4 S (u> = 0)/g 2 < x 2 > vs. temperature T at various 
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couplings. 

We see that as g 2 < x 2 > increases, the ratio falls 
sharply below the weak coupling value and indeed ulti- 
mately £ (u> = 0) becomes proportional to (g 2 < x 2 > 
) 1,/2 i-c. roughly to T 1 / 2 . The crossover may be un- 
derstood analytically from Eq |[ At low T and small 
g,a ~ it so one may expand in gx/a, generating the 
familiar Migdal series. At high T and large g a is 
negligible and the integral is dominated by a pole at 
x ~ it/g; taking account of the normalization of P(x) 

yields G/ oc ~ l/g(x 2 } 1 ^ 2 and so £ ~ glx 2 ) 1 ^ 2 . Thus 
the key to saturation is that at strong coupling the scat- 
tering rate continues to increase but at a rate less rapid 
than that given by second order perturbation theory. A 
similar result for the self energy of a model of carriers 
coupled to spin fluctuations was presented in ]lq] . 

We now add disorder scattering. Qualitatively one may 
think of impurities as adding an extra T-independent 
term to (x 2 ), so (x 2 ) -> (x 2 phonm {T) + x 2 mvunt ^ . At 
low T and small A the rates add, leading to Matthiessen's 
rule p = pphonon + Pimpurity At high T, physical 

1 /2 

quantities depend on (s 2 ) so the relative correc- 
tion due to impurities is of order x 2 mpunty /x 2 phonon {T): 
Matthiessen's rule does not apply. Fig 5 shows the re- 
sults of a direct calculation confirming this idea. A more 
sophisticated (non-CPA) treatment is required to repro- 
duce the observed negative dp/dT at very strong disor- 
der. 
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Fig 5 Impurity dependence of resistivity; A = 2.34, n=l/2. 

To summarize, we have shown that a simple model of 
electrons coupled to classical phonons and to static disor- 
der reproduces very well the essential features of the phe- 
nomenon usually described as resistivity saturation. The 
cause of the phenomenon was found to be a breakdown 
of the Migdal relation between lattice distortion ampli- 
tude and electron self energy. Comparison of our results 
for different carrier densities along with the band-mixing 
arguments of Ref indicates that the value of A/B will 
depend strongly on model details, so detailed compari- 



son of our calculation to experiments on specific materials 
is inappropriate, but material-specific calculations using 
the ideas and formalism put forward here would be of 
great interest. 

Our results imply that the term 'saturation' is a mis- 
nomer: there is no intrinsic maximum value of the high- 
temperature resistivity. Indeed at very high tempera- 
tures one expects classical diffusion with diffusion con- 
stant D vanishing or tending to a constant as T — > 00, im- 
plying p increasing indefinitely with T . This was found in 
'retracable path approximation' studies of the Hubbard 
model |L6| and recently in a model of carriers coupled 
to magnetically correlated spins |17| . From this perspec- 
tive the 'absence of saturation' |7[] in correlated electron 
materials is not by itself surprising; the interesting is- 
sue is the apparent smoothness of the crossover from the 
low-T coherent transport regime to the high-T classical 
diffusion regime. 
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